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Abstract. We investigate a relation between the Brownian motion on the unitary group and 
the most natural random walk on the symmetric group, based on Schur-Weyl duality. We use 
this relation to establish a convergent power series expansion for the expectation of a product of 
traces of powers of a random unitary matrix under the heat kernel measure. This expectation 
turns out to be the generating series of certain paths in the Cayley graph of the symmetric 
group. Using our expansion, we recover asymptotic results of Xu, Biane and Voiculescu. We 
give an interpretation of our main expansion in terms of random ramified coverings of a disk. 

1. Introduction 

In this paper, we are concerned with the asymptotics of large random unitary matrices dis- 
tributed according to the heat kernel measure. This problem has been studied first about ten 
years ago by P. Biane pQ and F. Xu [2]. It shares some similarities with the case of unitary 
matrices distributed under the Haar measure, studied by B. Collins and P. Sniady [3J 0]. The 
origin of our interest in this problem is the hypothetical existence of a large iV limit to the 
two-dimensional U(N) Yang-Mills theory. This limit has been investigated by physicists, in par- 
ticular by V. Kazakov and V. Kostov [5] and by D. Gross, in collaboration with W. Taylor [6], 
A. Matytsin [7] and R. Gopakumar [8]. In [9], I. Singer has given the name of "Master field" to 
this limit, which still has to be constructed. A. Sengupta has described in [10] the relationship 
between Yang-Mills theory and large unitary matrices. We refer the interested reader to this 
paper and will not develop this motivation further. Sengupta's work also contains some results 
whose study was at the origin of this paper (see Proposition 12.21 and the discussion thereafter). 

Our approach relies on the fact that the Schur-Weyl duality determines a (non-bijective) 
correspondence between conjugation-invariant objects on the unitary group on one hand and 
on the symmetric group on the other hand. To be specific, let n,N > 1 be integers. Let 
Pn : N '■ S n x U(N) — > GL((C N )® n ) be the classical representation. The set V n ,N of partitions 
of n with at most ./V parts indexes irreducible representations of both 6 n and U (N). If A is such 
a partition, let x X (resp. Xx) denote the corresponding character on S n (resp. U{N)). Let Z 
be an element of the centre of C[6 n ]. Let D be a conjugation- invariant distribution on U(N). 
Then the equalities 

(.1) VA G P n N , - — — — , 

X A (id) Xx{In) 

where X\(D) = Dxx, imply p n ^(Z <g> 1) = ^^(1 (g> D). The main observation, implicit in [6], 
is the following: the element Z = — J2i<k<l<n(^) ^ C[(5 n ] and the distribution D on 

U(N) defined by Dip = ^Ajjrmcp^N), where ^u(N) is the Laplace operator, satisfy ([I]). Now 
Z is, up to an additive constant, the generator of the most natural random walk on & n and it 
follows from this discussion that this random walk is closely related to the Brownian motion on 
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the unitary group. This relation is stated precisely and proved in Section 2. It is also partially 
generalized to the orthogonal and symplectic groups. 

In Section 3, we prove our main result, which is the following. 

Theorem 1.1 (see also Thm [3~3|) . Let N,n>l be integers. Let {Bt)t>o be a Brownian motion 
on U(N) starting at the identity and corresponding to the scalar product (X, Y) \— > —Ti(XY) on 
u(N). Let a be an element of & n . Let mi, . . . ,m r denote the lengths of the cycles of a. Then, 
for all t > 0, we have the following series expansion: 



For all T > 0, this expansion converges uniformly on (N,t) E N* x [0, T}. 

The coefficients S(a, k, d) count paths in the Cayley graph of the symmetric group & n . More 
specifically, we consider the Cayley graph of & n generated by all transpositions. For all ir € & n , 
we denote by \ir\ the graph distance between \tt\ and the identity. Then S(a, k, d) is the number 
of paths starting at a of length k and finishing at a point tt such that |7r| = |c| — (k — 2d). In 
particular, S(o~,k,d) = if \k — 2d\ > n : for each d > 0, the contribution of order N~ 2d is a 
polynomial in t. 

The coefficients S(a, k, d) depend on a only through its conjugacy class and can be expressed 
in terms of the representations of the symmetric group. In fact, Theorem 11.11 can be proved 
directly using the representation theory of the unitary and symmetric groups. We present this 
proof in Section 4. It is more systematical than the proof presented in Section 3 and should be 
easier to generalize, as also suggested by the work of Gross and Taylor [6]. 

The tools of representation theory allow us, in Section 5, to compute S(o~, k, d) when a is a 
cycle of length n. The expression involves Stirling numbers and it could hardly be called simple. 
Nevertheless, it allows us to count for all integer p the number of ways to write the cycle (1 . . . n) 
in & n as a product of p transpositions. 

In Section 6, we use our expansion to describe the asymptotic distribution of unitary matrices 
under the heat kernel measure as their size tends to infinity, thus recovering a result of P. Biane 
PQ. We also recover a result of F. Xu [2] on the asymptotic factorization of the expected values 
of products of traces. In order to describe the asymptotic distribution, we must compute the 
coefficients S(a, k, 0). The factorization result mentioned above reduces the problem to the case 
where a is an n-cycle. Unfortunately, the expression of ^((l . . . n), k, 0) obtained in Section 5 is 
not obviously equal to what it should be according to Biane's results. Thus, we compute this 
coefficient in a different way by using the relations between the geometry of the Cayley graph 
of the symmetric group and the lattice of non-crossing partitions. Then, in Section 7, we apply 
the same ideas related to non-crossing partitions and use Speicher's criterion of freeness to prove 
the asymptotic freeness of independent unitary matrices under the heat kernel measure. 

Finally, in Section 8, we give an interpretation of our formula in terms of random ramified 
coverings over a disk, thus proving a formula described by Gross and Taylor [6]. We define 
a probability measure on a certain set of ramified coverings over the disk and prove that the 
expectation computed in Theorem II .11 is the integral of a simple function - essentially N raised 
to a power equal to the Euler characteristic of the total space of the covering - against this 
measure. From this point of view, our expansion deserves to be called a genus expansion. 

It is a pleasure to thank Philippe Biane for several enlightening conversations. 
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2. Probabilistic aspects of Schur-Weyl duality 

In this first section, we establish formulae which relate the heat kernel measures on U(N), 
SU(N), SO(N) and Sp(N), to natural random walks in the symmetric group and the Brauer 
monoid. 

2.1. The unitary group. Let n and N be two positive integers. There is a natural action 
of each of the groups U(N) and & n on the vector space (C N )® n , defined as follows: for all 
U G U (N), a € & n and x\, . . . ,x n G C^, we set 

U ■ [x\ ® . . . ® x n ) = Ux\ ® . . . ® Ux n , 

(3) a ■ (xi ® . . . ® x n ) = x 0-1(1) ® . . . ® x a -i( n y 

It is a basic observation that these actions commute to each other. In particular, they determine 
an action p UtN of & n x U(N) on (C w ) 8n by 

Pn,N(cr, U)(xi ® . . . ® X n ) = Ux^-l^ ® . . . ® [7x CT -i (n) . 

Definition 2.1. Lei Mi, . . . , M n be N x N complex matrices. Let a be an element of & n . We 
denote by p u t (M\, . . . , M n ) the complex number 

p s J(Mi, ...,M n ) = Tr (C i V) ®n ((Mi ® . . . ® M n ) o PritN (a, I N )) 

= H Tr(M h ...M ir ). 

c=(h...i r ) 
cycle of a 

We set p s J{M) = p s J{M, . . . , M). 

The upper index st indicates that we use the standard trace rather than the normalized one 
in the definition. The letter p stand for "power sums", since p a t (M), as a symmetric function of 
the eigenvalues of M, is the product of power sums corresponding to the partition determined 
by a. Observe that, by definition, the character of the representation p n n is the function 

x Pn ,A^u)= Pa \u). 

The core result of Schur-Weyl duality is that the two subalgebras of End((C iV )® n ) generated 
respectively by the actions of U(N) and 6 n are each other's commutant. Let us explain why 
this makes a relation between the Brownian motion on U(N) and some element of the centre of 
the group algebra of <3 n unavoidable. 

Let u(N) denote the Lie algebra of U (N), which consists of the N xN anti-Hermitian complex 
matrices. Let U(u(N)) denote the enveloping algebra of u(N), which is canonically isomorphic 
to the algebra of left-invariant differential operators on U(N). Let also C[S n ] denote the group 
algebra of & n . The representation p n> N determines a homomorphism of associative algebras 
C[6 n ]®W(u(7V)) — ► End((C iV )® n ). The centre Z(u(N)) of W(u(7V)) is the space of bi-invariant 
differential operators on U(N). Since |O n ,jv(l ® Z(u(N))) commutes with p n ^(l,U) for every 
U G U(N), the Schur-Weyl duality asserts in particular that 

PniN {i ® z(u(N))) c PriiN (c[e n } ® 1). 

We are primarily interested in the Laplace operator, which is defined as follows. The IR-bilinear 
form (X,Y) = Tr(X*Y) = -Tr(XY) is a scalar product on u(N). Let (X u . . . , X N ?) be an 
orthonormal basis of u(N). Identifying the elements of u(N) with left-invariant vector fields 
on U(N), thus with first-order differential operators on U(N), the Laplace operator Amjv) is 

the differential operator YliLi X?- It corresponds to the Casimir element ® -^i °f the 

enveloping algebra of u(N). This element is central and does not depend on the choice of the 
orthonormal basis. Hence, ^u(N) is weu defined and bi-invariant. The discussion above shows 
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that, in the representation p n N, the Laplace operator of U(N) can be expressed as an element 
of C[(5 n ]. This is exactly what the main formula of this section does, in an explicit way. 

Let T n be the subset of & n consisting of all transpositions. We set A@ n = —^y^- + Yl T eT n T - 
The formula for the unitary group is the following. 



Proposition 2.2. For all integers, n,N > I, one has 

1 A \ Nn + n(n — 1) 
2 

Before we prove this formula, let us derive some of its consequences. 



(4) Pn,N (&<S n 1 + 1® ^U(N)) 



Proposition 2.3. For each a € & n , the function pff : U(N) — ► C satisfies the following 
relation: 

(5) \^U(N)P S a = ~^ ~ £ P. 



2 

T£T n 



st 

UT ' 



More generally, let Mi, . . . , M n be arbitrary NxN matrices. Then, regarding p^(MiU, . . . , M n U) 
as a function of U € U(N), one has 

(6) ±A u{N) p s J(MiU, M n U) = -^p*J(MiU, M n U) - £ p s J T (MiU, . . . , M n U). 

reT n 

Proof - Recall that p s J(M 1 U, M n U) = Tr ((Mi ® ... <g> M n ) o p n>N (cr, U)). Let us use the 
shorthand notation M = Mi ® . . . ® M n . We have 

lA u(N)P s J(M l U,...,M n U) = Tr(Mo PntN (a,U)op njN (l®~A u{N) )) 



Nn + n(n — 1) st 



2 

p s *(M x XJ, . . . , M n U) - Tr(M o p njN (a, U) o p„,7v(A en ® 1)). 



2 

The result follows immediately from the definition of A@ n . □ 

The function p s * depends only on the cycle structure of a. In concrete terms, if the lengths of 
the cycles of o are mi, . . . , m r , then p^(U) = Ti(U mi ) . . . Tr(U mr ). This redundant labelling is 
however nicely adapted to our problem, as equation ([5]) shows. Let us spell out the right hand 
side of this equality. The permutation a being fixed, the cycle structure of ar depends on the 
two points exchanged by the transposition r. If they belong to the same cycle of a, then this 
cycle is split into two cycles. A cycle of length m can be split into a cycle of length s and a 
cycle of length m — s by m distinct transpositions, unless m = 2s, in which case only ^ of these 
transpositions are distinct. If on the contrary the points exchanged by r belong to two distinct 
cycles of a, these two cycles are merged into a single cycle. Two cycles of lengths m and m! can 
be merged by mm' distinct permutations. Altogether, we find the following equation, which was 
already present in papers of Xu [2] and Sengupta [TO] . 



A u{N) (Tt(U mi )...Tr(U mr )) = -NnTi{U mi )...Ti{U r 



mi— l 



+ E"V rr (£ /mi ) • • -Tr([/™i) • • .Tr(£/ m '') ^ ^{U s )Ti{U m ^ s ) 

8=1 8=1 

r 

+ E ^i^jTr(^ mi )...Tr|L^0--- Tl {^O--- Tr (t /mr ) Tr(C/ m ^ 
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A remarkable feature of (J2J) is the fact that the element of C[S n ] which appears has coefficients 
of the same sign on the elements which are not the identity. Hence, up to an additive constant, 
it can be interpreted as the generator of a Markov chain on & n . This leads us to the following 
simple probabilistic interpretation of (j3]). 

Let us introduce the standard random walk on the Cayley graph of the symmetric group 
generated by the set of transpositions. It is the continuous-time Markov chain on © n with 
generator A@ n , that is, the chain which jumps at rate Q) from its current position a to err, 
where r is chosen uniformly at random among the Q) transpositions of & n . 

If a is a permutation, we denote by 1(a) the number of cycles of a. For example, r is a 
transposition if and only if £(t) = n — 1. 

Proposition 2.4. LetN,n > 1 be integers. Let (Bt)t>o be a Brownian motion on U(N) starting 
at the identity and corresponding to the scalar product (X, Y) i— > — Tr(XY) on u(N). Let (irt)t>o 
be a standard random walk on the Cayley graph of the symmetric group & n , independent of 



Nn + n(n—l) 



(Bt)t>o. Then the process (e 2 p**(_Bj)j is a martingale. In particular, 



w— 1 Nn-\-n{n — 1) , 

(7) E[rf t (Bt)]=e s 'E 



More generally, let M±, . . . , M n be arbitrary NxN complex matrices. Then the stochastic process 

( Nn + n(n-l) . \ 
e 2 p s ^ t (M\Bt, . . . , M n B t )j is a martingale and 

r j. -1 Nn-\-n(n — l) , r , 

(8) E [p% (M 1 B t , M n B t )} = e 2 * E [p s J Q (Mi, ... , M„)] . 

Proof — The process (tt, B) is a Markov process on & n x U(N) with generator A@ n (8 1 + 1 
\^U(N)- Consider the function p : & n x U(N) -> C defined by p(a, U) = p^(MiC7, . . . , M n U). 
By Proposition 12.3} this function satisfies the relation 

f . 1 A \ Nn + n(n - 1) 
( A Sn 1 + 1® - A U{N) J p = j '-p. 

( Nn+n(n-l) . \ 
e a X t {M 1 B t ,..., M n B t ) J is a martingale follows immediately. The 

last assertion follows from the fact that l?o = ^JV a-s. □ 



Let us turn to the proof of Proposition 12. 21 



Proof of Proposition 12.21 — The action of u(N) on (C )® n extends by complexification to 
gl(N, C) = u(JV) © iu(N). Let (X\, . . . , ^^2) be a real basis of u(N). It is also a complex basis 
of Ql(N, C). Define a N x N matrix g by 5^ = — Tr(XjX,). Since — Tr(- ■) is non-degenerate 
on gl N (C), the matrix g has an inverse the entries of which we denote by g l K Then it is 
easy to check that the element ^ . =1 g % ^Xi®Xj of the enveloping algebra is independent of the 
choice of the basis. Of course, by choosing our original basis of u(N) orthonormal, we find that 
this element is simply Amjn. 

In order to compute p n ,A r (l <8> Amjv)), we prefer to use another complex basis of gl(N, C) = 
Mjv(C), namely the canonical basis (^ij)ije{i jv}- F° r this basis, fl^jw = —SjkSu and 5 = g -1 . 
Hence, in the enveloping algebra of gl(iV, C), A (7 ( A r) = - 
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First, notice that /9 n> jv(l ® Eij)(x\ ® . . . ® x n ) = Y7k=i x i ® ■ • • ® Eij( x k) <8> • • • <8> x n . Hence, 

(AT \ AT 

1 ^ = 2 ^ ^ Id^" 1 ® ^ ® Id®'"*" 1 ® Eji ® Id ™-'- 1 + 

ij'=l y i>i=l KkKn 

AT n 

+ Id ^ 1 ® Id®" - * -1 . 

i,j=l k=l 

The last term is simply Nn times the identity. For the first part of the right hand side, observe 
that Ylij=i Eij ® Eji G End((C Ar )® 2 ) is the transposition operator x ® y i— ► y ® x, that is, the 
operator p2,AT ((12), ijv)- Finally, we have found that 

-Pn,Af(l ® ^U(N)) = Nn Id + Pn,N((kl),I N ). 

l<k^l<n 

The result follows. □ 

The results of this section still hold, after a minor modification, when U(N) is replaced by 
SU(N). Indeed, the orthogonal complement of su(N) in u(N) is the line generated by -4=ijy. 

Since /9 n ,Af(l ® (iiv ® -Tat)) = ^ 2 Id, the Casimir operator of su(iV) satisfies the relation 

Pn }N (l ® &SU(N)) = Pn,Ar(l ® A^)) + —Id. 

This modifies only the exponential factors in ([7]) and ([8]). 

We will explore further consequences of Proposition 12.41 in the rest of the paper. For the 
moment, we derive similar results for the orthogonal and symplectic group. 



2.2. The orthogonal group. Let us consider the action of SO(N) on (C*) 8 ™ defined by 
analogy with ([3]). The action of & n still commutes to that of SO(N), but, unless n = 1, 
the subalgebra of End((C )® n ) generated by the image of C[S n ] is strictly smaller than the 
commutant of the image of SO(N). Let us review briefly the operators which are classically 
used to describe this commutant. We denote by {ex, . . . , ejv} the canonical basis of C^. 

Definition 2.5. Let P be a partition of {1, . . . , 2n} into pairs. Define p n ,N{P) £ End((C^)® n ) 
by setting, for all . . . , i n G {1, . . . , N}, 

Pn,N{(3)(e h ® ...®e in ) = ^ ]J S hn e *n+i ® • • • ® e *2n- 

i n+ i,...,i 2 „e{i,...,Ar} {fc,«}e/3 

Observe that the partition {{1, n + 1}, . . . , {n, 2n}} is sent to the identity operator by p n .N- 
Let B n denote the set of partitions of {1, . . . , 2n} into pairs. The composition of the operators 
Pu,n{I3) corresponds to a monoid structure on B n which is easiest to understand on a picture. 
An element of B n is represented in a box with n dots on its top edge and n dots on its bottom 
edge. The dots on the top are labelled from 1 to n, from the left to the right. The dots on the 
bottom are labelled from n + 1 to 2n, from the left to the right too. A pairing is then simply 
represented by n chords which join the appropriate dots. Multiplication of pairings is done in 
the intuitive topological way by superposing boxes and, if necessary, removing the closed loops 
which have appeared. 
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Figure 1. Multiplication of two diagrams in the Brauer monoid. 

The monoid B n is called the Brauer monoid and its elements are called Brauer diagrams. 
The group & n is naturally a submonoic0 of B n , by the identification of a permutation a with 
the pairing {{1, <r(l) + n}, . . . , {n, cr(n) + n}}. The identification of & n with a subset of B n 
is compatible with our previous definition of Pn,jv m the sense that Pn,N{&) is the same if we 
consider a as a permutation or as a Brauer diagram. 

The correct statement of Schur-Weyl duality in the present context is that the subalgebras 
of End((C w )® n ) generated by SO(N) and B n are each other's commutant (see Let Pn,N 

denote the morphism of monoids 

p n)N : B n x SO(N) — GL({C N )® n ). 

Just as in the unitary case, this action determines a morphism of associative algebras p n N '■ 
C[B n ]®U($o(N)) — ► End((C Ar )® n ). 

By analogy to the unitary case, let us define "power sums" functions associated to Brauer 
diagrams. Given (3 £ B n and Mi, . . . , M n G Mjv(C), set 

pf (Mi, . . . ,M„) = Tr((Mi ® . . . ® M n ) o ^(/J)). 

In particular, the character of p n> Ar is given by Xp„ jv(A ^) = P s a{^)- 

The number pi* (Mi, . . . , M n ) is a product of traces of words in the matrices Mi,*Mi , . . . , 
M n ,*M n . Let us describe in more detail how to compute p s ^(In)- Let (3 be a Brauer diagram. 
Consider the graph with vertices {1, . . . , n} and unoriented edges {k, I}, where k and I are such 
that there exist k' £ {/c, A; + n} and /' 6 {/, / + n} with {k' , /'} G /3. This is the graph obtained 
by identifying the top edge with the bottom edge in the graphical representation of (3. Then 
each vertex has degree 2 in this graph. Hence, it is a union of disjoint unoriented cycles. If [3 
belongs to & n C B n , this cycle structure is of course that of f3 as a permutation, apart from the 
orientation which is lost. In general, let £(f3) denote the number of cycles in this graph. Then 
p%(I N ) = N^\ 

Let us define an element of C[-B n ] as follows. Given k and I two integers such that 1 < k < 
I < n, we define the element (kl) of B n as the following pairing: 

(kl) = {{k, I}, {n + k,n + I}} U |J {{i,n + i}}. 

ie{l,...,n}-{k,l} 



In fact, S„ C B n is exactly the subset of invertible elements. Indeed, for [3 € B n , let T((3) be the set of pairs 
{k, 1} £ (3 such that 1 < k, I < n. In words, T(/3) is the set of chords in the diagram of /3 which join two dots on 
the top edge of the box. It is clear that T(/3i/32) 3 T(f3i) for all /3i,/?2 £ B n . Hence, T(/3) must be empty for (5 
to be invertible. More generally, it is not difficult to check that, given (3 and /3' in B n , there exists /3" £ B n such 
that 0(3" = f3' if and only if T{(3) C T((3'). 
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Figure 2. The elements (24) and (24) of B 6 . 

Let C n be the subset of B n consisting of all the element of the form (kl). We now define 
As n = h a. Thanks to the inclusion & n C B n , we still see Ag n as an element 

of C[-B n ]. The formula for the orthogonal group is the following. 
Proposition 2.6. For all integers, n, N > 1, one has 

(9) Pn,N (A Sn 1 + 1® A SO (jv)) = -- — + Pn,iV (A Bn ® 1) • 

Proof — The computation is very similar to that we made in the unitary case. Endow so(N) 
with the scalar product (X,Y) = —Tr(XY). The basis (^4y)i<i<i<jVj with Ay = — Eji, is 
orthogonal and (Aij,Aij) = 2 for all i < j. Hence, A SO (n) = \ Yli<i<j<N ^ij <8> Aij. We have 

p n , N (l ® A so(N) ) = E Id 0fe - 1 ®^®Id^- fc - 1 ®A ij ®Id 0n -' + 

l<fe<Kn l<i<j<N 



2 o j^jigin— fc 



+^E E h^- 1 ®^ 

fc=l l<i<j<N 
N 

Id® fc_1 ® £^ <8> Id®'"^ 1 ® £?y ® Id ™"' 

KkKn M=l 

- E E Id ^~ 1 ® ^ ® Id '~ fc ~ 1 ® ^ Id "~' - 2 Id 

l<fc<Kn i,j=l 



^ Pn,Ar(((fcO " (kl)) ® 1) - ^ Id. 



l<fc<Kn 

The result follows. □ 



The following proposition is proved just as Proposition 127 
Proposition 2.7. For a// /? G B n , the following relation holds: 

(10) A so{N)P % = -^^ P f - E * + E Pgr 

tST„ aeC„ 

More generally, let M\, . . . , M n be arbitrary NxN matrices. Then, regarding p^(MiR, . . . , M n R) 
as a function of R £ SO(N), 

(at — iw 

&SO{N)Pp(MiR, M n R) = - V - > pf(M l R, M n R) 

(11) - E Pfr(MiR, M n R) + Y, PU M i R > • • • > M nR)- 

reT n a&C„ 
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It seems more difficult to find a probabilistic interpretation of ([9]) than in the unitary case, 
because the element of C[-B n ] which appears does not have coefficients of the same sign on all 
elements not equal to 1. 

2.3. The symplectic group. Nothing really new is needed to treat the case of the symplectic 
group. Let us describe briefly the results. 

Let J € M2jv(C) denote the matrix ( ® n r ^he symplectic group is defined by 



-In 

Sp(N) = {S <E U(2N) : *5J5 = J}. It acts naturally on ((C 2JV ) 0n ). The action of the Brauer 
monoid needs to be slightly modified to fit the symplectic case. If (3 belongs to B n , we define 
the operator p n ,2N(P) by setting, for allii, {1, ... , 2iV}, 

Pn,2Jv(/3) {e h ® . . . <g> e in ) = ^ ]J J lkk e in+1 <g> . . . <g> e i2n . 

in+l,— ,»2ne{l, —,n} {k,l}Gf3 

Then we have an action p n ^N '■ B n x Sp(N) — > End((C 2Ar )® n ) and the images of B n and Sp(N) 
generate two algebras which are each other's commutant. 

The Lie algebra sp(N) is endowed with the scalar product {X, Y) = —Ti(XY) and we denote 
by Agp(jv) the corresponding Laplace operator. The main formula is the following. 

Proposition 2.8. For all integers, n,N > 1, one has 

(12) p n ,2N (A Sn ® 1 + 1 ® 2A Sp(N) ) = -(2N + l)n + p n , 2N (A Bn ® 1) . 

Proof — Just as in the unitary case, it is more convenient to use complexification. The Lie 
algebra sp(N, C) = sp(N) © isp(N) is the Lie subalgebra of g[(2N,C) defined by the relation 

l XJ = —JX. It consists of the matrices f ^ ^ , where A is an arbitrary N x N matrix 

and B, C are two symmetric N x N matrices. We use the following basis of sp(N, C): 

Mj = Eij - Ej + N,i+N 1 < i,j < N 

Bij = Eij +N + E j>i+ N 1 < i < j < N 

Cij = E i+N j + E j+N)i 1 < i < j < N 

Di = E u + N 1 < i < N 



D i+N = E i+Nti 1 < i < N. 

The bilinear form (■, •) takes the following values on this basis: 

(A ij ,A ji ) = -2 l<i,j<N 
(By,Ci i > = -2 l<i<j<iV 
(D u D t+N ) = -l l<i<N. 
The other values are zero. It follows that the Casimir element of sp(N, C) is equal to 

A sp(N,c) = ^2 Ai i ® A ^ ~\ ^2 ( Bi i ® Ci i + Ci i ® 

±<i,j<N 1<*<J<JV 

- ^2 (Di ® D i+N + D i+N ®Di). 

l<i<N 

The formula follows now by a direct computation. In order to recognize operators of the form 
(kl) and (kl), observe that, when n = 2 for example, 

2N 

P2,2n((12),I n ) = Eij^Eji, 
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N 

P2,2Jv((12),ijv) = ^ (Eij E i+ NJ+N + Ej+Nj+N <S> Ejj 
i,j = l 

—Eij+N ® E i+ Nj - E i+ N,j <S> Eij+x) ■ 

□ 

Proposition 2.9. For a// (3 £ B n , the following relation holds: 

(13) 2A 5p(iV) ^ = -(2JV + l)n j# - £ + £ pf a . 

reT n a&Cn 

More generally, let Mi, . . . , M n be arbitrary 2Nx2N matrices. Then, regarding p s p(M\S, M n S) 
as a function of R E Sp(N), one has 

2A 5p(7V) ^ i t (M 1 S,...,M n S) = -(2N + l)npf(M 1 S,...,M n S) 

(14) - ^^ t T (M 1 5,...,M n 5)+ £ p^(M 1 S,...,M n S). 

r&T n aeC n 

3. The power series expansion 

Let us denote by Tttv = 77 Tr the normalized trace on Mjv(C). Let Mi, . . . ,M n be N x N 
matrices. Let a be an element of & n . We denote by p a (Mi, . . . , M n ) the number 

p a (M 1 ,...,M n )= 11 Tr N (M h ...M ir ). 

c=(ii...i r ) 
cycle of cr 

We denote by 1(a) the number of cycles of a, so that p a = N^ e ^p^. 

In this section, we exploit the result of Proposition 12,41 and derive a convergent power series 

expansion of E p a (B±_) when B is a Brownian motion on U(N). This expansion involves 
L w J 

combinatorial coefficients, which count paths in the Cayley graph of & n . We start by discussing 
these paths and introducing some notation. 



3.1. The Cayley graph of the symmetric group. Fix n > 1. The Cayley graph of & n 
generated by T n can be described as follows: the vertices of this graph are the elements of & n 
and two permutations a\ and 02 are joined by an edge if and only if aia^ 1 is a transposition. 
It is a fundamental observation that, if a% and o~2 are joined by an edge, then £(a\) and t(a%) 
differ exactly by 1. Indeed, multiplying a permutation by a transposition splits a cycle into 
two shorter cycles if the points exchanged by the transposition belong originally to the same 
cycle, and otherwise combines together the two cycles which contain the points exchanged by 
the transposition. 

A finite sequence (o"o, . . . , Ofc) of permutations such that a% is joined to <7j + i by an edge for 
each i € {0, . . . , k — 1} is called a path of length k. The distance between two permutations is 
the smallest length of a path which joins them. This distance can be computed explicitly as 
follows. 

Let us introduce the notation \a\ = n — 1(a). We have \a\ € {0, . . . , n — 1} and \a\ = (resp. 
1, resp. n — 1) if and only if a is the identity (resp. a transposition, resp. a n-cycle). Other 
values of |<r| do not characterize uniquely the conjugacy class of a. It is well-known and easy to 
check that |o"| is the minimal number of transpositions required to write a. In other words, the 
graph distance between two permutations a\ and 02 in the Cayley graph is given by | c"^ - 1 c"2 1 - 
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It turns out that the paths which play the most important role in our problem are those 
which tend to get closer to the identity. Let 7 = (<ro, . . . be a path. Recall that, for all 
i G {0, . . . ,k — 1}, one has £(<Ti+i) = t{pi) db 1. We call defect of 7 and denote by (£(7) the 
number of steps which increase the distance to the identity. In symbols, 

d{i) = #{ie{O,...,fc-l}:|0- <+ i| = |<7 < | + l} 

= #{i G {0, . . . , k - 1} : £(a i+1 ) = £(a { ) - 1}. 

The following lemma is straightforward. 

Lemma 3.1. Let 7 = (<7q, . . . , a^) be a path. Then 2d(*y) = k — (£(a k ) ~ £( a o))- 

For a, a' in & n and k > 0, let us denote by II^ (cj — ► a') the set of paths of length k which 
start at a and finish at a' . Let us also denote by Tiki®) the set of all paths of length k starting 
at a and by Il(cr — > a') the set of all paths from a to a'. Notice that the cardinality of 11^(0") is 
equal to (9) ■ Let us finally define the coefficients which appear in the expansion. 

Definition 3.2. Consider a G & n and two integers k,d>0. We set 

S(a,k,d) = #{7 eIL k (o) : d( 7 ) = d}. 

In words, S(a,k, d) is the number of paths in the Cayley graph of 6 n starting at a, of length k 
and with defect d. 

Observe that the adjoint action of & n on itself determines an action of & n on its Cayley graph 
by automorphisms. Thus, S(a, k, d) depends only on the conjugacy class of a. 

3.2. The main expansion. 

Theorem 3.3. Let N,n>l be integers. Let (Bt)t>o be a Brownian motion on U(N) starting 
at the identity and corresponding to the scalar product (X, Y) 1— > —Tr(XY) on u(N). Let 
Mi, . . . ,M n be arbitrary N x N complex matrices. Let a be an element of & n . Then, for all 
t > 0, we have the following series expansions: 
(15) 

E[p a (M 1 Bt_,...,M n Bt_)]=e-f J2 2 #Tl k {o ^ a') p a ,{M u . . . ,M n ). 

k,d=0 ' |cr'| = |cr|-A:+2d 

In particular, if mi, . . . ,m r denote the lengths of the cycles of a, then 

(16) E[Tr N (Bp...Tr N (BT)] = e ~x £ ^-J^ S (a,k,d). 

X X k,d=0 

For all T > 0, both expansions converge uniformly on (N,t) € N* x [0, T]. 

In order to understand the role of the defect of a path in our problem, let us write down the 
result corresponding to Proposition 12.31 for the functions p a . As explained earlier, the number 
of cycles of ar can be either 1(a) + 1 or 1(a) — 1, respectively when the two points exchanged 
by r belong to the same cycle of a or to two distinct cycles. For each permutation a G & n , we 
are led to partition T n into two classes F(a) and C(a), those which fragment a cycle of a and 
those which coagulate two cycles. More precisely, 

F(a) = {r G T(n) : i(ar) = 1(a) + 1} and C(a) = {r G T(n) : £(ar) = £(a) - 1}. 

The following result is now a straightforward consequence of Proposition [ 
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Proposition 3.4. Let a be a permutation in & n . Let Mi, . . . , M n be N x N matrices. Then 
the following relation holds: 
1 



—AuwpAMx U,..., M n U) = -- 



p a (M 1 U,...,M n U) 



+ P*AM 1 U,...,M n U) + ^ E Por(M 1 ?7,...,M n C0 



reF(cr) 



According to this result, each step which increases the distance to the identity is penalized 
by a weight iV -2 . In the proof of the power series expansion, we use the following lemma. 

Lemma 3.5. Let t > and N > be real numbers. For all a, a' € & n and e 6 {—1, 1}, define 

-k+k 



k=0 



Then the matrices (Mj: . ; ) (T(7 /ge n and (M o .^,) - ]Cr / e g n are eac/i other's inverse. 
Proof — Let us define an endomorphism L of C[(S n ] by setting, for all / € C[S n ], 

{Lf)(a)= Y, /M + 4f E /^ r )- 



reF( CT ) 



reC((7) 



One checks easily that the matrix M^ a , is the matrix of the operator e on C[S n ] and the 
result follows. □ 



Proof of Theorem 13.31 Consider T > 0. We claim that the right-hand side of (I15p is a nor- 
mally convergent series on (N, t) € N* x [0, T\. Indeed, let us define K = max{|p (7 (Mi, . . . , M n )\ : 
a E 6 n }. Then, for all N > 1 and all t € [0, T], the sum of the absolute values of the terms of 
the series is smaller than 



+OQ rpk +2° 



Ke-^^2 — Y J S(a,k,d)=Ke 



n(n — 2) , 



k 

k=0 ' d=0 

The assertion on the uniform convergence of the expansions follows. 

In order to prove (|15p . we start from the expression given by Proposition 12.41 a t time and 
with an arbitrary deterministic initial condition ttq = a. It reads 



(17) 



pi\ (MiB±, . . . ,M n B± 



7T = a 



nt _ n(n—l)t 



e 2 



—p*{M u ...,M n ). 



We expand the left hand side by using the properties of (^t)t>o- This chain jumps at rate and 
its jump chain is a standard discrete-time random walk on the Cayley graph of & n , independent 
of the jump times. Thus, the left-hand side of (fT7]) is equal to 



E 

fc=0 



_(n\±_ I n 

e \2> n 



k\N k ( n \ k ^— ' ^— ' 



(n\K t 
\2) <J' 



p s a t ,(M 1 B ± ,...,M n B ± ) 



ee„ 7 en fe (o-^cr') 



where the expectation is now only with respect to the Brownian motion. After simplification 
and switching to normalized traces, (|17p becomes 



00 +k 



Vex € & n , ^ E \p°> (MtBt_ M n Bt_ )] X) jfei 



t k #U k (a -> a') 



fc=0 



e 2 P(T (Afi, 



,M n 
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We recognize the expression of M^ a , and, by Lemma E3J we conclude that for all a £ 6 n , 



E 



p a (MiBt, M n Bt) 



e 2 



fc=0 



A:! 



The first formula follows from the fact that \a\ 
yields the second formula. 



n - ((a). Setting Mi, 



, M^v equal to /at 
□ 



Most of the coefficients which appear in the expansion (|16p are zero. More precisely, the 
situation is the following. 

Lemma 3.6. Let 7 be a path of length k and defect d starting at a. Then the following inequal- 
ities hold: 

0<d<k and 2d - {1(a) - 1) < k < 2d + (n - ((a)). 
In particular, \k — 2d\ < n — 1. 

Moreover, let d > be given. Then S(a, 2d + (n — 1(a)), d) > and, if d > 1(a) — 1, then 
S{a, 2d - (£(a) -l),d)>0. Finally, if d < 1(a) - I, then S(a, d, d) > 0. 

Proof — Assume that the path finishes at a k . Then the first two inequalities reflect simply the 
fact that 1 < i{a k ) < n. 

To prove the second part of the statement, consider d > 0. Recall that a is fixed. Let us 
construct a longest possible path starting at a with defect d. For this, we minimize the defect 
at each step. First, we build a path by going from a down to the identity through a geodesic. 
This takes n — £(a) steps and the defect of the path is still zero. Then the path must make one 
step up. Immediately after this, it can go down to the identity again. It can repeat this at most 
d times without its defect becoming larger than d. By then it has length 2d + (n — t(a)). Thus 
we have constructed a path of length 2d + (n — 1(a)) with defect d. A similar argument works 
for a shortest path of given defect. □ 

In particular, for all d > 0, the contribution of order N~ 2d to ¥,[p a (B±)] is a polynomial 

function of t of degree 2d + (n — ((a)) and in which the smallest exponent of t is max(d, 2d — 
(((a)-l)). 

3.3. Examples. Let us work out explicitly a few examples. 

For n = 1: there is a single path in the Cayley graph of ©1. It has length and defect 0. Thus, 
we recover the well-known formula 



E 



Tr N (B 



e 2 . 



For n = 2: for each k > there is a unique path of length k starting at the identity. It has 
defect L^J- Thus, 

00 (—i) k t k ( t 1 t N 

= e~* — ■ — -rrrr- = e~* f cosh — — — sinh — 



E 



Tr N (B±) 

N 



k\N 2 ^i 



k=0 



N N 



N 



Similarly, for each k > 0, there is a unique path of length k starting at (12). It has defect [_f 
Thus, 



E 



Tr N (B\) 



E 

k=0 



k+k 



( cosh 



N 



N sinh ■ 



N 



For n = 3: the situation is a bit more complicated but it is still possible to compute everything by 
hand. For a path starting at the identity of length k and defect d, we must have 2d — 2 < k < 2d. 
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Hence, if k is odd, it must be equal to 2d — 1. So, for all I > 1, S^id, 21 — 1,1) = 3 . If is 
even, then two situations are possible. We leave it as an exercise to check that, for all / > 1, 
S(id,2l,l) = 3 21 ' 1 and 5(id,2/,/ + 1) = 2.3 2l -\ Finally, 5(id,0,0) = 1. We find 



E 



Tr N (B 



1 + 



JV 2 + 2 
3JV 2 



cosh 



3i 

JV 



1 . , 3t 

N Smh N 



Similarly, we find 



E 



Tr N (B\)Tr N (B_ 



2 



cosh 



3i 
JV 



N 2 + 2 3* 

— — sinh — 

3JV JV 



E 



1MB 



3t 
2 



1 + 



N 2 + 2 



cosh 



3i 

JV 



1 



JV sinh 



3i 

JV 



For n > 4, it seems difficult to determine all the coefficients at once and by hand. Nevertheless, 
the following diagram, which indicates how many edges join the various conjugacy classes of 64 
in the Cayley graph allows one to compute specific values of S(a, k, d). 




Figure 3. The Cayley graph of S4 modulo conjugation. 

For instance, one can use it to prove the following formulae : 
Tr N (B ± ) 4 



e 2t E 



1 + ^ (-6t + 3t 2 ) + ^(15i 2 - 20i 3 + 5i 4 ) + O(-^g), 



e 2t E 



Tr N (B\) = (1 - 6i + 8t 2 

N 



h + ^ (10t2 



58 



71 



-t° + —r 

3 4 



16 c; N 1 , 

T * 5 )+0(^). 



3.4. The case of SU(N). Let us conclude this section by stating without proof the following 
analogue of Theorem 13.31 in the case of the the special unitary group. This theorem is proved 
exactly like its unitary version, by using the observation made at the end of Section 12.11 

Theorem 3.7. Let JV, n > 1 be integers. Let (Bt)t>o be a Brownian motion on SU(N) starting 
at the identity and corresponding to the scalar product (X,Y) 1— > —Tr(XY) on su(JV). Let 
M\, . . . ,M n be arbitrary JV x JV complex matrices. Let a be an element of & n . Then, for all 
t > 0, we have the following series expansion: 

+00 



E 



Pa{M t B 



± ,...,M n B ±/ 

N N 



nt_ 1 n t 

e 2 2n' 2 



E 



-1)H 



k-ik 



k\N 2d 



(18) 

In particular, if mi , 



k,d=0 

#Il k (a^o-')p a ,(Mi, 



E 

\a'\ = \a\-k+2d 

, m r denote the lengths of the cycles of a, then 



Mr, 



(19) 



2 +00 (_-i\kfk 
E[Tr N (Bp...Tr N (Bp] = e-^ ^ 

k,d=0 



S(a, k, d). 
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For all T > 0, both expansions converge uniformly on (N,t) G N* x [0, T]. 



4. A REPRESENTATION-THEORETIC DERIVATION OF THE POWER SERIES EXPANSIONS 

In this section, we give an alternative derivation of the expansions (|16p and (I19p . based on the 
representation theory of the unitary and symmetric groups and the relations between symmetric 
functions. This approach is less elementary than the one adopted in the previous sections but 
we believe that it is more likely to allow generalizations. In Section O we will use it to compute 
some of the coefficients S(a, k, d). 

4.1. Expansion for the unitary group. The integers N, n > 1 are fixed throughout this 
section. We write A h n if A = (Ai > . . . > A r > 0) is a partition of n. The integer r is called 
the length of A and we denote it by ^(A). We denote the set of all partitions by V. 

Let A be a partition of n. We denote by s\ the Schur function associated to the partition 
A, whose definition is given in |12t 1.3]. For all U 6 U(N), the number s\(U) is defined as the 
value of s\ on the eigenvalues of U. We will use the fact that, if £(X) > N, then the symmetric 
polynomial in N variables determined by s\ is the zero polynomial. This follows for example 
from the expression of s\ as a determinant in the elementary symmetric functions \12\ 1.3, (3.5)]. 

Recall the definition of the power sums, that is, the functions p s J: : U (N) — > C for a G & n 
(see Definition 12. ip . 

The Schur functions and the power sums are related as follows. Let x A : ®n — * C denote the 
character of the irreducible representation of & n associated with A. Then one has the following 
pair of relations [121 1-7,(7.7)]: 

(20) VA h n, s x = - £ X X M, 

o-eS„ 



(21) V,£6 n , p s J = Y,x\o-)s x . 

Ahn 

The set U(N) of isomorphism classes of irreducible representations (irreps) of U(N) is in 
one-to-one correspondence with the set of non- increasing sequences a = (or > . . . > a^v) of 
elements of Z. Even when some of the a^s are negative, the Schur function s a is well-defined 
and the character of the irrep a is Xa(U) = s a (U). 

Let (Bt)t>o be the Brownian motion on U(N) of Theorem 13.31 Let dU denote the normalized 
Haar measure on U(N). For each t > 0, let Qt denote the heat kernel at time t on U(N), that 
is, the density of the distribution of Bt with respect to the Haar measure. Our main result is 
the following reformulation of (|16p . 

Theorem 4.1. Let N,n>l be integers. Let a be an element of & n . Then, for all t > 0, 

f +oc (—^ k t k 

(22) N-^ / pZ*(E0Q$(E0 dU = e-T £ 

Ju i N ) k,d=0 

In the course of the proof, we admit two lemmas which we prove afterwards. We have pre- 
ferred this order to the strict logical order to make the proof easier to follow. 



Proof — If t = 0, the result is clearly true. When t > 0, the proof consists in expanding Qt into 
the sum of its Fourier series and turning all quantities related to U(N) into quantities related to 
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& n . For all t > 0, the function Qt is smooth on U(N) and invariant by conjugation. It admits 
the following uniformly convergent Fourier expansion |13L Thm 4.4]: 

(23) Q t (U) = Y, e"^s«(ljv)^), 

aeZ L 

where the number 02(a) is characterized by the equality Au(N)Xa = — C2(a)Xa- Using the 
relation (|2ip to expand p^(U), we find the following expression for the left-hand side of ([22]): 

f P s J(U)Q ± (U)dU = V e-^SailNWfr) [ 7JU)s^U)dU. 



By the orthogonality properties of the characters of irreps, the integral in the right-hand side is 
zero unless a = fi. Hence, we can replace the sum over a and n by a sum over the partitions \i 
such that fj,\- n and £(n) < iV: 

(24) / p s J(U)Q,(U)dU = Yl e-^s^x^a). 

Ju W M hn,£ (At )<7V 

We still need to express s„(Jjv) an d C2(//) in terms of quantities related to the symmetric group. 

In order to compute s^(ijv), we use the relation Qgft. Let us define Q = E CT eS„ N e ^a. This 
notation is borrowed from [6]. Then (I20p implies the equality 

(25) S „(I*) = 

In Lemma 14.31 we will prove that X**(fi) = if £(p) > N. This allows us to drop the restriction 
£(n) < iV in the summation. 

Let us compute 02(11), the eigenvalue of A-u(N) associated to s^. Thanks to ([5]), we know 
the value of Au(N)Pa f° r an an< ^ p0|^ expresses s M as a linear combination of power sums. 
Combining these two equations, we find 

(26) *uw a v = - Nns » "^E E X^r)p s J, 

where T n is the set of the transpositions of <3 n . We now use the following consequence of Schur's 
lemma: whenever x belongs to the group algebra C[S n ] and y to the centre of the group algebra, 

x' M (x)x^(y) 



(27) V M hn, x^xy)- v;!jl . 

This relation implies that the last term of (|26p is equal to n(n — 1)— ^^y^s^. Hence, 



Y M (fl2)) 

(28) C2 ( M ) = niV + n(n-l)^-^. 

Combining $M$), and (J2SD, we find 

Jow * b k ' k n! ^ 2N xH1 > > 

By Lemma 14.21 below, the sum over [i is equal to N e ^ J2d>o N~ 2d S(a, k, d). The result follows 
immediately. □ 
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Lemma 4.2. Let a be an element of & n . Let k > be an integer. Then 

1 > ~ 2^ n \ I 2 J ■ 

d>0 fihn V 

Proof — There is no issue of convergence, since the sum on the left-hand side is finite. Let 
7 = (do, . . . , cifc) be a path of defect d. By Lemma l3~Tj £(cFk) = £{&$) + k — 2d. Hence, 

d>o o-'ee„ 

Now, #rifc((T — > o"') is the number of fe-tuples (n, . . . , r/%) G such that <rri . . .t^ = a 1 . A 
standard computation based on the fact that X^hn X /1 ( i )x A '( cr ) = "-^o-.id and on (p7|) leads to 

„ n , f\ ST X"(l)x"(g- V) ^(n - 1) x "((12)) \ k 

The result follows by summing over a' and applying (|27[) again. □ 



Lemma 4.3. Define SI G C[S n ] by SI = Ylaee„ ^ <?"■ Then the following relations hold. 

1. For all fi\- n such that < N, x M (^) = n!s M (ijv)- 

2. For all fi\~n such that £(//) > iV, x M (^) = 0. 

Proof — The first assertion follows immediately from (|20|) . 

In order to prove the second assertion, let us introduce the Jucys-Murphy elements X±, . . . , X n 
of C[6 n ], defined by Jfi = and X { = (1 i) + (2 i) + . . . + (i - 1 i) for i € {2, ... , n}. They 
generate a maximal Abelian subalgebra of C[S n ]. In particular, they can be simultaneously 
diagonalized in every irreducible representation of & n . We borrow the following statements 
from [14j . 

Let fx = (/ii, . . . ,fj, r ) be a partition of n. The subset = G (N*) 2 : j < fii} of Z 2 is 

called the diagram of [i. An element (i,j) of -D M is called a box and its content is defined as the 
integer c(i,j) = j — i. 

The space of the irreducible representation of 6 n associated to /x admits a basis which diago- 
nalizes the Jucys-Murphy elements and is indexed by the bijections t : — > {1, . . . , n} which 
are increasing in each variable. These bijections are usually called tableaux. The eigenvalue of 
the Jucys-Murphy element X^ on the vector associated to the tableau t is the content of the 
box t' 1 ^). 

We need also the following well-known fact: for every k G {0, . . . , n — 1}, the fc-th elementary 
symmetric function of the Jucys-Murphy elements is equal to J]o-e6„ •"•M=& a > where \a\ = 
n — i{a). This can be proved as follows. For each m G {1, . . . , n}, let us imbed & m into <S n as 
the subgroup which leaves {m+1, . . . , re} invariant. For all m G {1, . . . , re} and k G {0, . . . , m— 1}, 
set S mifc = Eaee™, J-kNfc °~- Let us also define e m , fc = Yli 1< ...<i k < m x h ■■■ x i k - We need to 
prove that e m ^ = S m This is clearly true if k G {0, 1}. The general case follows by induction 
on rre, each inductive step being proved by induction on k, thanks to the relations 

^"m,k = ^m~l,k + ^m—\,k—X X k an d &m,k = &m—l,k + 6m— l,fe— l X k- 

From the equality proved in the last paragraph and the relation \a\ = n — £(o~), we deduce 
the following equality in the polynomial ring C[© n ][z]: 

n 

(30) H(z + X t )= ^2 z e{a) °- 

1 = 1 crSSn 
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Evaluating at z = N and applying x^i we find 

(n \ n lii 

iiiN+Xi) = £ n( iv + c (*" i «))=x M (i)nn( iv +i-^ 
1=1 / t tableau i=l i=l j=l 

If > 2V+1, then (iV+1, 1) is a box of £> M , whose content is —TV. It follows that x M (^) = 
in this case. □ 



4.2. Expansion for the special unitary group. Let us apply a similar analysis to the special 
unitary group in order to derive (|19p in another way. 

By restriction, any irrep of U(N) determines an irrep of SU(N) and the restrictions of a, a' G 
Z^ are isomorphic if and only if there exists k € Z such that a' = {ot\ + k, . . . , ajv + k). Hence, 
the set of irreps of SU(N) is in one-to-one correspondence with the set of partitions of length 
at most N — 1 and the character of the irreducible representation corresponding to a partition 
A is given by the Schur function s\. 

Let (Bt)t>o be the Brownian motion on SU(N) of Theorem [19j Let dU denote the Haar 
measure on SU(N). For each t > 0, let Qt denote the heat kernel at time t on SU(N), that is, 
the density of the distribution of B t with respect to the Haar measure. We reformulate (fT9|) as 
follows. 

Theorem 4.4. Let N,n>l be integers. Let a be an element of & n . Then, for all t > 0, 

(31) N-*W / rf*(l7)Q^(E0 dU = e-^+i^ £ i^-^M- 

J SU(N) k,d=0 ' 

Proof — If t = 0, both sides are equal to 1. Assume that i > 0. The Fourier expansion of Qt is 
then the following: 

(32) Qt(U)= Y, e-^sxiL^T^U), 

xev 
e(\)<N-i 

where now c' 2 (A) is defined by the equality A su ^ N ^xx = ~~ c^Wxx- Combined with the relation 
(|2TD . it implies 

/ P s J(U)Q ± (U)dU= V e-^sxilN^io-) [ JxJUjs^U) dU. 

JSU(N) N X( . Vlihn JSU(N) 

£{X)<N-1 

When £(fi) > N, is identically zero on SU(N). When £(fi) < N — 1, then the integral in the 
right-hand side of the last equation is equal to S\ !fl . Let us consider the terms of the sum for 
which l(pL) = N. In this case, let us write fj, = // + (/Xjv, . . . , /j>n), so that £(fj-') < N — l. Then the 
integral is equal to 5\^i and we may assume that A = //. In this case, s\{Ln) = s^i(Lj^) = s^Ln) 
and c' 2 (A) = c' 2 (fj-')- Hence, 

/ p s J(U)Q<,(U)dU = Yl e-^s,(L N ) X ^a)+ £ e-^ s,{L N ) X ^(a). 

Jsu ( N ) n\-n,l{^)<N-l n\-n,e(n)=N 

Let us first compute c 2 (/x) when £(n) < N — 1. For this, we use the fact that is an eigenvector 
of A[/(jv) whose restriction to SU(N) is s^. Since A SU ( N }p^ = (Ajj^-j + j?) p& whenever 
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n , we find, thanks to (1271) . 



&SU(N)S, = (-Nn- n{n - l)^rf + ¥ V 



, X^((12)) n 2 

When i(jx) = JV, we are interested in s^i but is still an eigenvector of Au(N) whose restriction 
to SU(N) is s^i. Hence, 

A A ( AT I n ^!((lg)) | " 2> \ / \ 

^S£/(JV)V = A>su(N)Si* = [~Nn - n(n - 1) ,^ + — 1 s M = -c 2 (/x )v 
Thus, in both sums, the argument of the exponential is — ^ + — t n ^ N ^ X ^(i)^ • Using this 



fact and the first assertion of Lemma 14.31 we find 

nt „2 t n(n-l) x M ((12)) X /' i ( ) : x I' I rr \ 



/ p s J(U)Q ± (U) dU = e~~ + ^ V e" 



'sum-- < ^m<n nl 

The second assertion of Lemma 14.31 tells us that we can remove the restriction £(/j,) < N, since 
the other terms are zero. After expanding the exponential, Lemma 14.21 allows us to finish the 
proof just as in the unitary case. □ 



5. Computation of 5((1 . . .n),k, d). 



In this section, we apply the methods of representation theory to the computation of some 
of the coefficients which appear in our main expansions, namely the coefficients S((l . . . re), k, d) 
for all n,k,d > 0. 

Let us recall the definition of the Stirling cycle numbers, or Stirling numbers of the first kind 
s(n, k), also denoted by [?] . They are characterized by the identities in C[x] 



x(x — 1) . . . (x — n + 1) = 
valid for all re > 0. In other words, 



k=0 



x 



(-l) n - k e n . k (l, 



(-1) 



n—k 



E « 

l<il<...<i n _fc<n— 1 



■ 1<n—ki 



where e n _£ denotes the (n — fe)-th elementary symmetric function. 

By applying the alternating character to the identity (|30p . we find the relation 



e{cj)x 1 ^ = x(x 



l)...(x-n + l) 



E 

fe=0 



from which we deduce that | [^] | is the number of elements of S n with exactly k cycles, or 
in other words at distance n — k from the identity. In particular, [q] = 0. Let us make the 
convention that [?1 = if k < 0. The main result of this section is the following. 



Proposition 5.1. For all n,k,d> 0, 

S((l...n),k,d) = - ~ 
re 

r,s,l,m>0 
r+s=n— 1 
l+m=n— 1— k+2d 



\l+r 



r!s! 



n 

2 (S 



r 


s + 1 




r + 1 




s + 1 - l_ 




r + 1 — m 
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Proof — Instead of computing S({1 . . . n), k, d) we compute the sum of S(a, k, d) when a spans 
the set of all n-cycles. Dividing the result by (n — 1)! yields S((l . . . n), k, d). Now, a path of 
length k starting at an n-cycle has defect d if and only if it ends at a distance n — 1 — k + 2d from 
the identity. Let us recall some of the notation used in the proof of Lemma 14.31 The integer n 
being fixed, we set E r = ^|o-|=r a - Hence, 

S((l...n),k,d) = j^-^ Yl S(cr,k,d) 

£reSn,|c|=n— 1 



(33) 



(n 


-1)1 




1 


(n 


-1)! 




1 


(n 


-1)! 



Y #n fe (cT - TT) 

CT,7r£S n ,|f |=n— 1, 1 7r | =ri — 1— k+2d 



E 

Ahn 



n! V X A (id) 



Now we use the following fact, which is a consequence of the description of the representations 
of & n given by Okounkov and Vershik [T3] and recalled briefly in the proof of Lemma [ 



(34) 



Vr e {0, ...,n- 1} , 



X A (id) 



e r ({<:(□) : □ S A}) . 



In words, the right-hand side of this equation is the r-th elementary symmetric function of the 
contents of the boxes of the diagram of A. In particular, if the diagram of A has at least two 
boxes of content 0, then x A (X n _i) = 0. Hence, the non-zero terms of the sum (j33|) arise from 
the partitions which are hooks, that is, of the form rj r = (n — r l r ) for some r G {0, . . . , n — 1}. 
This fact is well known (see for example the appendix of [E]), and rj r is the representation 
/\ r St, where St is the restriction of the natural representation of & n on C n to the hyperplane of 
equation {z\ + . . . + z n = 0}. This representation is of degree [ n ~ ) and x Vr ((^ ■ ■ • n )) = ( — l) r - 
Let us introduce the notation s = n — 1 — r. It follows easily from (|34p that 



(35) 



Vr G {0, ...,n- 1} , 



X^(Si) _ n(n- 1) 
X^(id) ~ 2 



nr 



n 

2 {S 



The contents of the boxes of r\ r are {— r, . . . , 0, . . . , s}. Hence, by the definition of the Stirling 
numbers and (|34p . we have for all r £ {0, . . . , n — 1} 



X Vr (id) 



e n _i_ fc+2( i(-'r, ... ,0, ... ,s) 

^ e/ (l,..., S )(-ire m (l,...,r) 



-m=n— 1— fc+2d 



(36) 



E 



«,m>0 
Z+m=n— 1— k+2d 



s + l 




r + 1 


s + l-l_ 




r + 1 — m 



Finally, combining (j35[) and ()36p . we find the expected result. 



□ 



It seems that Proposition 15.11 should allow one to find a simple generating function for the 
numbers ^((l . . . n), k, d). Our best result in this direction is the following. We use the notation 
(x) n = x(x — 1) . . . (x — n + 1). 
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Proposition 5.2. For all n,N >0, one has 

E(-l) k t k , , (-1)™ 1 -1- (s-r) 2 n 2 t 2 (fit . \ (fit, \ 

LJ^SCP... „),*,,() = <->- YJ e^< » ( T (.-r)J ( - T (-r)j 

k,d>0 r,s>0 v / s+1 \ / r 

r+s=n—l 

We emphasize that this generating function is, unfortunately, exponential with respect to d 
instead of k. 

In Proposition 15.11 when k takes the largest possible value given n and d, namely n — 1 + 2d, 
then I and m must be equal to and the identity r?l = 1 simplifies greatly the expression. This 
leads us to the following corollary. 

Corollary 5.3. Let n > 1 be an integer. For each p > 0, Zei c njJ) denote the number of distinct 
ways in which the cycle (l...n) € & n can be written as a product of p transpositions. The 
number c n . p is non-zero if and only if p = n — 1 + 2d for some d > 0. In this case, 

c n , p = S((l...n),n-l + 2d,d) = 

For eac/i n > 1, one /ias iZie equality 

Cn,P^ = -e 2 x (l-e nx ) . 
p>0 ^ 

In particular, c n) „,_i = n n_2 , c^+i = ^(n 2 — l)n n+1 and 

Cn,„+ 3 = 7^(5n - 7)(n + 3)(n + 2)(n 2 - l)n n+3 . 
57oU 

Remark 5.4. T/ie waZne of c n ^ n -\ is classical. The sequence (c n . n +i) n >i zs known as A060603 
in the Online Encyclopedia of Integer Sequences [16] . 




6. Asymptotic distribution 

One of the consequences of Theorem 13.31 is that the limit as N tends to infinity of E[p ff (Sj.)] 
exists. Using Lemma 13.61 we get the following result. 

Proposition 6.1. Consider a € & n . The limit of IL[p a (B±)] as N tends to infinity exists and 
it is equal to 

lim eLb,)1 =e -fV(-lf^. 

k=0 

Unfortunately, Proposition l5.1l does not seem to lead easily to a simple expression for S(a, k, 0) 
nor even S((l . . . n), k, 0). In this section, we determine a simple expression of S(a, k, 0) for all 
a and k > 0. For this, we prove a factorization property and use the relation between the 
metric geometry of the Cayley graph of & n and the lattice of non-crossing partitions of the cycle 
(1 . . . n). The fact that the two expressions of S((l . . . n), k, 0) given by Propositions 15. l l and 16.61 
agree is not obvious, at least for the author. 
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6.1. The factorization property. The factorization property is the following result. It re- 
duces the problem of the determination of S(a, k, 0) to the case where a is a cycle. 



Proposition 6.2. Let mi, . . . ,m r be positive integers. Then 



lim E 

AT— >oo 



Tv N (BT)...Tr N (BT) 

N N 



lim E 

JV— >oo 



N 



lim E 



More precisely, 



(37) 



E 



Tr N (B 



mi \ 

_L / 

JV 



E 



Tr N (BT] 



...E 



Tt n (BT 

N 



Tr N (B[ 



0(N~ 



uniformly in t on bounded intervals. 



We start by proving the following property of the numbers S(o~, k, 0). It is in fact equivalent 
to the proposition. 

Proposition 6.3. Consider a E & n . Assume that a = c\. . . ci( a ) is the decomposition of a as 
a product of cycles with disjoint support. Then 



(38) 



V/c>0, S((r, k, 0) 



kl 



h+...+i 



III... l e ^\ 



S(ci,h,0) . . . S( c *(er)>^(<T)>0). 



Proof — The number S(a, k, 0) is the number of paths of length k starting at a and which at 
each step move towards a permutation with one more cycle than their current position. As we 
already observed several times, each step of such a path corresponds to the multiplication by a 
transposition which exchanges two points which belong to the same cycle of a. There is thus 
a natural partition of the set of all steps of such a path, according to the cycle of a in which 
their support is contained. Let us introduce some notation. Let (ctq = a,a\, . . . ,ak) be a path 
with defect zero. For each i E {1, . . . , k}, set Tj = o~~\oi. Let {C\, . . . , Q(cr)) be the partition of 
{1, . . . ,k} determined by the fact that i E Cj if and only if the support of Tj is contained in the 
support of Cj. Then it is clear that, for all j € {1, . . . ,£(a)}, the transpositions {ji,i £ Cj) are 
the steps of a path with defect zero starting from Cj. 

Hence, constructing a path of length k starting at a and with defect zero is equivalent to 
constructing a collection of £(a) paths with defect zero starting at c\ , . . . , cir a \ respectively, 
whose lengths h,. ■ ■ , ^( CT ) add up to k, and a shuffling of the steps of these paths, that is, a 
sequence (Ci, . . . , CV^)) of subsets of {1, . . . , k} which partition {1, . . . , k} and whose cardinals 
are h,... ,le( a ) respectively. The equation (|38|) is just the translation in symbols of the last 
sentence. □ 



Proof of Proposition 16.21 — By Theorem 13.31 and Proposition 16.31 tli e terms of degree 
of the difference on the left hand side of (I37p vanish. Hence, this difference is of the form 
N~ 2 F(t, N~ 2 ) for some entire function F. The result follows. □ 

We have observed after Definition 13.21 that S(o~,k,d) depends only on the conjugacy class of 
a. Hence, we need to compute S((l . . . m), k, 0). The arguments of the proof of Proposition 
16.31 show that the paths of defect starting at (1 . . . m) stay in & m if we identify & m with the 
subgroup of & n which leaves {m + 1, . . . , n} invariant. Hence, we are reduced to the computation 
of S((l . . . n), k, 0) for all n > 1 and k E {0, . . . , n — 1}. This computations involves non-crossing 
partitions. For the sake of being self-contained, we give a brief review of the properties of 
non-crossing partitions that we use. 
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6.2. Non-crossing partitions. Let P = {Pi, . . . ,Pf\ be a partition of {1,. . . , n}. The par- 
tition P is said to be non-crossing if there does not exist i,j,k,l G {l,...,re} such that 
i < j < k < I and r, s G {1,... ,£} with r ^ s such that i, k G P r and j,l G P s . Another 
way to formulate the fact that P is non-crossing is the following. For each class Pj of the parti- 
tion, let Hj denote the convex hull in C of {e « : k G Pj}. Then P is non-crossing if and only 
if for all i, j G {1, ...,£}, i ^ j => Hi n Hj = 0. This notion is relative to the cyclic order on 
{1, . . . ,n} determined by (1 . . . n). We denote by NC(n) the set of non-crossing partitions of 
the cycle (1 . . . n). This set has been first considered by Kreweras in |17j . 

The fineness relation between partitions restricted to NC(n) makes NC(n) a poset. More 
precisely, we say that P\ =^ P2 if every class of P\ is contained in a class of P<i- The poset 
(NC(n), =4) is in fact a lattice, which means that suprema and infima exist. There is in particular 
a maximum, {{1, . . . ,n}}, which we denote by l n , and a minimum, {{1}, . . . , {fi}}, which we 
denote by n . The poset NC[n) can be made into a graph by joining two partitions P and Q if 
they are distinct and comparable, say P -< Q, and the interval [P, Q] = {R G NC(n) : P ^ R =^ 
Q} is reduced to {P, Q}. 

A non-crossing partition P = (Pi, . . . , Pi) of the cycle (1 . . . n) determines an element op of 
S n as follows: take the cycles of op to be the classes of P with the cyclic order induced by 
(1 . . . n). In symbols, if i G Pj, then 

ap(i) = (1 . . . n) k (i) , where k = min{/ > 1 : (1 . . . n) l (i) G Pj}. 

In particular, <7o n = id and a± n = (I...71). The partial order on NC(n) corresponds via the 
mapping P 1— > op to the following partial order on ©,„. 

Consider (J\,oi G S n . Recall that |<7i| = n — l(o\), the minimal number of terms of a 
decomposition of o\ in a product of transpositions, is the distance from id to o\ in the Cayley 
graph of & n generated by T n . By definition, we say that o\ =^ 02 if |cr 2 1 = + o"2|- In 
words, o"i ^ 02 if and only if there exists a geodesic path from id to 02 through o\ . The identity is 
the minimum of & n for this partial order, and the n-cycles the (pairwise incomparable) maximal 
elements. The next lemma is well known and its proof is left to the reader. 

Lemma 6.4. The mapping from NC(n) to & n which sends a partition P to the permutation 
ap is an isomorphism of posets from NC(n) onto [id, (1 . . . n)\ = {a G & n : a ^ (1 . . . n)}. 

As a consequence of this Lemma, ^((l . . . n), k, 0) is the number of decreasing paths of length 
k starting at l n in NC(n). It turns out to be easier to count increasing paths in NC(n) starting 
at n . They are in one-to-one correspondence by the following duality property of NC(n) 
discovered by Kreweras. 

For a G [id, (1 . . . n)], let us introduce K(a) = . . . n). It is readily checked that K is a 

decreasing bijection of [id, (1 . . . n)\. The corresponding decreasing bijection of NC(n) is called 
the Kreweras complementation. It can be described combinatorially at the level of non-crossing 
partitions as follows. 

Given a partition R of {1, . . . ,n} and a partition S of {1, . . . , n} ~ {1, . . . , n}, let R U S 
denote the partition of {1, 1,2,2, ... ,n, n} obtained by merging R and S. Even if R and S are 
non-crossing, R U S may be crossing with respect to the cyclic order (1, 1,2,2, ... ,n, n). Now 
let P be a non-crossing partition of {1, . . . ,n}. The partition K(P) is by definition the largest 
element of NC(n) such that P U K(P) is non-crossing. 

The following result summarizes this discussion of non-crossing partitions in relation to our 
problem. 

Proposition 6.5. For all n > 1 and k > 0, S((l . . . n), k, 0) is the number of increasing paths 
of length k starting at {{!}, . . . , {n}} in the lattice of non-crossing partitions of (1 ... n). 
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Figure 4. The Kreweras complement of {{1, 3, 12}, {2}, {4, 8, 9}, {5, 6, 7}, {10, 11}} 
is {{1, 2}, {3, 9, 11}, {4, 7}, {5}, {6}, {8}, {10}, {12}}. 

It remains to count these paths. To do this, we use the fact that 5((1 . . . n),n — 1,0) = n n ~ 2 . 
This is a classical result of combinatorics, since ^((l . . .n),n — 1,0) is the number of ways to 
write a n-cycle as a product of n — 1 transpositions. It is also a special case of Corollary 15.31 

Proposition 6.6. For all n > 1 and k > 0, 



S((l...n),fc,0)=^ 1 jn*- 1 . 
It is understood that this number is zero if k > n. 

Proof — We count the increasing paths of length k in NC(n) starting at {{1}, . . . , {n}} by first 
regrouping them according to their terminal point. The possible terminal points of these paths 
are exactly the non-crossing partitions of (1 . . . n) into n — k classes. Such partitions may be 
classified according to the number of singletons they contain, the number of pairs, and so on. 

Let si, . . . , s n be non-negative integers such that si+. . .+s n = n—k and s\+2s2+- ■ -+ns n = n. 
We say that a partition is of type (si, . . . , s n ) if it contains exactly Sj classes of cardinal i for each 
i G {1, . . . , n}. The number of non-crossing partitions of type (s±, . . . , s n ) has been computed 

by Kreweras 1 171. It is equal to — — ■ -. 

(k + l)!si! . . . s n \ 

Let P be a non-crossing partition of type (si, . . . , s n ). An increasing path from {{1}, . . . , {n}} 
to P looked at in the reverse direction is a decreasing path from P to {{1}, . . . , {fi}}- There 
are as many such paths as there are geodesic paths from the permutation ap induced by P to 
the identity, that is, S(ap,k, 0). Let us apply Lemma [6731 to compute this number. The only 
non-zero term in the sum corresponds to the situation where l{ = |cj| for each i £ {1, . . . , £(o~p)}. 
Since S(c, |c|,0) = m m ~ 2 for every cycle c of size m, we find the formula 

S((l...n),k,0)= y M - t 2° S2 3 1S3 ...n("- 2 ^, 

u h ' ; ^ (k + l)!si! ...s n \ l\ S2 . . . (n — l)! Sn 

(s 1 ,...,s n ) 

where the sum is extended to all possible types of partitions of {1, . . . , n} into n — k classes. Let 
us enumerate the possible types by enumerating the partitions themselves. If we do this, each 
type will appear as many times as the number of partitions of this specific type. The number 
of partitions of type (si, . . . , s n ) is 



l! si si! . . . n\ Sn s n \ 
Hence, we find the following expression: 



5((l...n),fc,0) = ^ T ^l( 1 ^2 



(2-1)52 _ _ n (n~l)s 
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where the sum runs over all partitions of {1, . . . , n} into n — k classes. Now the right hand side 
has been computed by Kreweras in [18] and it is equal to 

s( ( i... n ), M ) = -A-f n ~ 1 ) n k = ( n \ n k-i 

u k + l\n-k-l) \k + y 

This is the expected result. □ 

Let us state separately the following result which has been used in the course of this proof. 
Lemma 6.7. Let P G NC(n) be a partition of type (s\, . . . , s n ). There are exactly 

(s 2 + ... + (n-l)g n )! 20S2 n (n-2)s n 
l! s 2 . . . (n - 

increasing paths from n to P. 

Remark 6.8. The explicit expression of the large N limit of the moments of B_t_ obtained in 

N 

this section has been stated by Singer in [9] and proved by Biane in pQ . 

The asymptotic distribution of B± as N tends to infinity is the unique probability measure 

N 

fit on the group U of complex numbers of modulus 1 such that, for all n G N ; J v z n fi(dz) = 

f v z~ n fit(dz) = e~~ Ylk=o ^~nk\ (fc+i) • This expression of the moments of fit is not very easy 
to handle, if only numerically, because it is an alternated sum of large numbers. Let us point 
out two analytical ways of studying fit- 

The S -transform of fit is its moment generating function, defined by St(z) = f v dfi t (£) = 
Y^n=i m n,tz n , where m n; t is the n-th moment of fit- Since \m ri: t\ < 1 for all n and t, the function 
St is holomorphic on the unit disk D = {z : \z\ < 1}. Moreover, since St(z) = z + 0(z 2 ), there 
exist a reciprocal function to St in a neighbourhood of 0, which we denote by Xt- It turns out 
that xt is much simpler than St: Xt(z) = ^rj-e^ z+5 \ as one can easily check by using Lagrange's 
inversion formula. To put it more concisely, the measure fit is fully characterized by the following 
relation, valid for z in a neighbourhood o/0: 

/ 7—^- d ^(0 = ! + *■ 

iul- -^je tz e2^ 

By studying xt, Biane proved in jT5] the following facts. For each t > 0, the measure fit has a 
continuous density with respect to the uniform measure on U. For t G (0,4], this density is zero 
exactly on the set 



jo 



e 



\9 — 7r| < 2arctan J V*(4 



t 1 



t 2 



For t > 4, the density of fit is positive. Finally, for all t > 0, the density of fit at e ld is a real 
analytic function of 6 on the relative interior of its support. 

Another way of studying fit is to observe that m n j = \ irm § e~ ntz (l + ~) dz, the integral 
being along any contour of index 1 with respect to 0. In |, Gross and Matytsin use this 
expression and the saddle point method to exhibit the following phase transition with respect to 

3 

t: if t G (0,4), m n> t decays with n like n~2 , whereas for t G (4, +00), it decays exponentially. 

_ 4 

Along the same lines, one can check that for t = 4, m n ^ decays like n 3 . This indicates that the 
density of fi^ is less regular than the density of fit fort G (0,4). This behaviour is consistent at 
a heuristical level with a general result of Biane about additive convolution with the semi-circle 
law [20]. 
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6.3. Almost sure convergence. The material gathered so far allows us to prove very easily 
the following result. 

Proposition 6.9. Consider a £ & n . Then, uniformly in t on bounded intervals, 



Var 



0(N~ 



In particular, on any probability space on which a Brownian motion on U(N) is defined for N 
large enough, the following convergence holds almost surely and in L 2 : 



lim Pa (B±) = e-^ V (-l)*S(a,k,0y-. 



k=0 



Proof — Let us denote by a x a the element of &2n which sends i on a(i) and n + i on n + o~(i) 
for each i € {1, . . . , re}. With this notation, p 2 = p^xa, so that 

2 



Var [p*(Bt_) 

{-l) k t k 
k\N 2d 



E 



E 



E 

k,d=o 



S(a x a, k, d) 



kl 



S(a,h,di)S(a, h,d 2 ) 



h+l2=k,di+d,2=d 

where the last expression follows after simplification from Theorem 13.31 In the term correspond- 
ing to d = 0, both d\ and d 2 must be equal to 0. By the same argument of support as in the 
proof of Proposition 16.31 we find 

S(axa,k,0)= nn S (^h,0)S(a,l 2 ,0). 



h+h=k 



h\l 2 \ 



The result follows immediately. 



□ 



7. Asymptotic freeness 

In this section, we prove that independent Brownian motions on U(N) converge in distribu- 
tion, as N tends to infinity, towards free non-commutative random variables. We do not consider 
*-freeness, that is, we do not consider products involving B~ t . In fact, the asymptotic *-freeness 

of independent Brownian motions follows from a general result of Voiculescu (see [21] and [1] 
Lemma 6] for details). Here we use Speicher's characterization of freeness by the vanishing of 
mixed free cumulants. This combinatorial characterization is very well suited to the approach 
we have adopted in this paper. 

7.1. The factorization property. Let us start by slightly improving Proposition 16.21 by in- 
cluding the extra deterministic matrices M\, . . . , M n of Theorem [331 We consider these matrices 
as elements of the non-commutative probability space (Mn(C), Tr^r) and speak of their distri- 
bution accordingly. 

Proposition 7.1. Let (m[ N \ . . . , Mn)^>i be a sequence of families of N x N matrices. 
Assume that this sequence converges in distribution. Consider a € & n . Write a as a product of 
cycles : a = (i l>x . . . i x>mi ) . . . {it(a),\ ■ ■ ■ ^( CT ), m£(<r) ) ■ Then 

(39) lim E P Jm\ N) B t M^ N) B t ) = TT lim E Tr N ( M^ N) B t ... M} N) B 
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More precisely, 

(40) E[p a (Mi N) B j _,...,M^B j _)] -]jETr N ^B^...M^lB j _)=0(N-% 

r=l 

uniformly in t on bounded intervals. It is understood that all limits exist. 
Proof — According to Theorem 13.31 the left hand side of (|39H is equal to 

(41) ^E^f £ # n ^--0 ^pAm[ n \...M n) )- 

k=0 \cr'\ = \cr\-k 

This last limit exists for all a' by the assumption that the family (Afj , • • • , Mn ) converges in 
distribution. All permutations a' which contribute to the sum satisfy on one hand lu'l = \o~\ — k, 
hence Ic/cr -1 ] > k, and #rifc(<7 — > a') > 0, hence Icr'cr - 1 1 < k. Hence, only permutations a' such 
that a' =^ a contribute and (HTT) can be rewritten as 



(42) e-f £ ( \ n kW| (aW) lim p ff (Mf», . . . , M™). 

Let ci, . . . , Q( CT ) denote the cycles of c It is not difficult to check that the interval [id, a] in the 
poset S n is isomorphic to the product of intervals Y\^=l IM, c r] by the mapping (ai, . . . , a^^) i— > 
ai . . . a^). Consider a' ^ a and write <r' = «i . . . accordingly. Then |cjV _1 | = laic^ 1 ) + 
... + |«^(<t)% ct ) I- Moreover, by the same argument of shuffling used in the computation of 
S((l...n),k,0), 



I / -in *M 

(7 (7 I' 



#n kV -!|(a -» a') = — H #n |arC -i!( Cr - a r ). 



H r v = i |a r c r |! ,, : , 

Hence, (I42D can be written as 



(43) r[ e -^F^ y ^ V # n, -!.(*-> a,) lim Tr w (Mf } . . . M \ N) 

r=l a r =?c r 1 1 



and this is just the right-hand side of (|39p . 

The equation (|40p follows from (|39p just as in Proposition 16.21 □ 



7.2. Free cumulants. As a preliminary to the proof of the asymptotic freeness, we compute the 
free cumulants of the limiting distribution of B± . Let (A, (p) be a non-commutative probability 
space and ut an element of A such that B± converges in distribution to ut as N tends to infinity. 

N 

We have spent a substantial part of this paper proving that the moments of ut are given by 




Given a G A and a € 6„ with cycle lengths (mi, . . . ,m^r a \), let us use the notation tp a {a) = 
(p(a mi ) . . . ip(a me< -< T )). The free cumulants of ut form a family of complex numbers (k n (ut)) ne \j >i © 
and they are characterized by the identity 

(45) Vn > l,Vcr g & n , <fcr{ut) = k a '(ut). 

<t'=<<7 



28 THIERRY LEVY 

Let us use the notation k n = kn mfl y It is an elementary property of the free cumulants that 
they are multiplicative, in that k a = k mi . . . k mt{<y) when (mi, . . . , m f ^) are the cycle lengths of 
a. 

Proposition 7.2. The free cumulants of ut are given by 

, , v nt (-nt)™" 1 
46 k n (u t ) = e~ \ \ 

n\n — 1)\ 

More generally, if a € & n , then 

(47) k a (u t ) = e~f iz^#n H (id^a). 

|CJ|! 

Proof — Let us put <p(uf) under the form of the right hand side of (I45p . Applying Theorem 
13.31 using Kreweras complementation and using Lemma 16.71 we h n d 

nt _ f_^k(i...t») _1 l 

¥>(«?) = e~ E U n) - 1 |, #n k (i...n)-i|((l---n)^^) 
(74(1.. .n) 1 V ' " ' '' 

<T=^(l...n) 

(e-T(-t)Ol-i)^i(e-T (-t)^ ) 8 * . . . (e~f (-t)"- 1 ^- 2 )^ 

0' S1 1 |S2 fn — IV s ™ 
cr^(l.. .n) v 

where si, . . . , s n are respectively the number of fixed points of a, and the numbers of transposi- 
tions, 3-cycles, . . ., n-cycles in the decomposition of a. By comparing this expression with (|45f) . 
we find the desired expression for the cumulants of ut- □ 

Let us recall briefly Speicher's characterization of freeness by the vanishing of mixed free 
cumulants |22j . Let a±, . . . ,a n be non-commutative random variables on a space (A, (p), where 
tp is a tracial state. For a € & n , the number m CT (ai, . . . , a n ) is defined by 



m CT (ai, . . . ,a n ) = \ [ ^( a ii • • • a v)- 



c cycle of a 
c=(h...i r ) 

It is well defined thanks to the fact that <p is tracial. The numbers m CT (ai, . . . , a n ) are the mixed 
moments of a%, . . . , a n . The relation 



m CT (ai, ... ,a n ) = ^ K'(ai, ... ,a n ) 



a> Re- 



characterizes the family of numbers k a (ai, . . . , a n ), the mixed free cumulants of a±, . . . ,a n . 

Speicher's characterization of freeness is the following. Let (dk)k>i be a family of elements of 
A. Then this family is free if and only if, for all n > 2 and all i\, . . . , i n > 1 such that i r ^ i s 
for some r,sE {1, . . . ,n}, fc(i... n ) (a^ , . . . ,a in ) = 0. 

Theorem 7.3. Let (B^ N ' k ^)^^>i be a family of Brownian motions, such that, for all N, k > 1, 
£>(N,k) 2 S a Brownian motion on U(N) and, for all N > 1, the Brownian motions (B^ N ' k ^)k>i 
are independent. Let (£&)&>! be a sequence of non-negative real numbers. 



SCHUR-WEYL DUALITY AND THE HEAT KERNEL MEASURE ON THE UNITARY GROUP 29 
Then, as N tends to infinity, the family of non- commutative random variables (B^' k ^)k>i 

N 

converges in distribution towards a free family (b^)k>i of non-commutative random variables 
such that, for all k > 1, b^ k ' has the distribution of Ut h . 

Proof — We prove the result for finite families {B^ N>i,ke{i,—,K} f° r some finite K, by 
induction on K. The case K = 1 is settled by our computation of the asymptotic distribution 
of Sj , that is, Propositions 16, l l and 16.61 

Let K > 2 be an integer and let us assume that the property is proved for K — 1 independent 
Brownian motions. We need to prove that the mixed free cumulants of B^' 1 ^, . . . ,B^J K ' tend 

JV N 

(N k) 

to zero as N tends to infinity. We regard B\ k ' as elements of the non-commutative probability 

w 

space (L°°(Q, P)®Mjv(C), E(g>TrAr) and we use the notation m a and k a accordingly. In particular, 
with our previous notation, m a =Ep (7 . 

What we need to prove is that, for all n > 2, all a G & n , all i±, . . . ,i n € {1, . . . , K} not all 
equal, 



lim m ff (Bp l) ,...,Bp" ) )= V lim k a ,{B^\ . . . , b[ N '^ 

Vre{l,...,n}, i a , {r) =i r 



By the factorization property (j39[) and Fubini's theorem, the left hand side is multiplicative with 
respect to the cycle decomposition of a. The right hand side is also clearly multiplicative, hence, 
it suffices to consider the case where a = (1 . . . n). In this case, we are looking at the expected 

trace of a product B^' 11 ^ . . . B[ N ' in ^ where at least two factors are distinct. Of course, the case 

TV" Tv" 

where one of the K possible factors does not appear is treated by the induction hypothesis. Let 
us assume that the K factors appear, in particular B^- N ' l \ Up to cyclic permutation, which does 
not affect its trace, the product above can be put under the form W\ . . . Wr Bt^ , 

AT W 

for some r > 1 and some products W[ , . . . , Wy of factors among B^ N,2 \ . . . , B^ N,K \ Our 
previous results show that 



lim ETr N (b[ N >^...B^A = e"f V 



CT <l...r) 



#n k(1 ... n) -i|((l . . . r) - <r) lim m^wf , . . . , W r (7V) ) 
J2 ^f^#n w (id^a) lim m* (ff) (wf°,...,wW) 

' — ' \<7 ! AT— >oo v ' 
cr^(l...r) 



V k a {uti ) hm m^wf , . . . , 

g ■ A/ — »nn 



N— >oo 
<x=$(l...r) 

where we have changed a in K(a) = . . . r) between the first and the second line. 

The term linijv^oo m K(a){w[ N \ ... , Wr N ^) is equal to a sum of limits of free cumulants of 

the factors appearing in W± N \ . . . , Wr N ^ in this order. By induction, only pure free cumulants 
appear, those which do not involve more than one Brownian motion in each cycle of the per- 
mutation. Moreover, by the combinatorial description of K(a) as a non-crossing partition, only 
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such cumulants appear that remain non-crossing when they are merged with a. In symbols, 



lim E Tr w ( Bp . . . B^A = V lim k„{Bf> h \ . . . , Bp n) ). 

V ^ w 7 <r<l...n) W n 

Vfee{l,...,n}, ifc=v (fc) 

This is exactly what we expected. □ 



8. Large N Yang-Mills theory on a disk and branching covers 

In this section, we explain how our main expansion relates the Brownian motion on the 
unitary group to a natural model of random branching covers on this disk. In doing this, we 
give a rigorous proof of results which are stated in [6] . 

Let D be the closed disk of radius 1 centred at the origin O of M 2 . Let n > 1 be an integer. 
Let A be a partition of n. We define the set lZ n> \ as the set of isomorphism classes of ramified 
coverings it : R — ► D which satisfy the following properties. 

1. R is a ramified covering of degree n. 

2. For each ramification point x € D of R which is not the origin O, R has a generic ramification 
type at x, in that #7r~ 1 (x) = n — 1. 

3. The monodromy of R along the boundary of D belongs to the conjugacy class of & n corre- 
sponding to A. 

An element of 7Z n .\ is allowed to be ramified over O, with any kind of ramification. The set of 
its ramification points distinct from O is called its locus of generic ramification. It is contained 
in the interior of D — {0}, which we denote by D* . 

Let X be a finite subset of D*. We define !Z n ,\,x as the subset of 1Z n ,\ formed by the coverings 
whose locus of generic ramification is X. 

The set 1Z n ,\,x is in natural one-to-one correspondence with a set of equivalence classes of 
paths in the Cayley graph of & n as follows. Assume that X = {x\, . . . ,Xk}- Choose a point b 
on the boundary of D. By the interior of a simple closed continuous curve based at b, we mean 
the bounded connected component of the complement of its range. Let C,C\, . . . , C& be simple 
closed curves in D based at b with pairwise disjoint interiors such that the interior of C contains 
O and, for each % G {1, . . . , k}, the interior of Cj contains Xj. We assume that this is done in such 
a way that the curve CC^ . . . C\ is homotopic to the boundary of D in D — (X U {O}). Then the 
monodromies ao, t\, . . . , of an element R € TZ n ,\,x along the curves C, C\, . . . , are defined 
in & n up to simultaneous conjugation and their orbit characterizes R. The assumptions made 
on R imply that n, . . . , t/. are transpositions and oot\ ■ ■ ■ Tp- belongs to A. 

Let V n \k be the set of paths in the Cayley graph of & n which start at an element of the con- 
jugacy class determined by A. The symmetric group acts on V n! \,k by conjugation. The mapping 
from TZ n> x,x which associates to R the orbit of the path (o"o r fc • • • T i, &o T k ■ ■ ■ T"2, ■ ■ ■ > &oTk, &o) 
is a bijection. Moreover, the cardinal of the stabilizer of this orbit is equal to the cardinal of 
the automorphism group Aut(i?) of R. Hence, the image on TZ nj \,x of the counting measure on 
V n ,\.k by the mapping V n xk — > V n ^^j& n ~ TZ n:X ,x is the measure 

PnXX = £ # Aut(i?) SR - 

This measure is finite and satisfies p n ,x,x(^) — Co) ■ 

There is a natural topology on lZ n \ , which is generated by the sets 



0(R,U) — {R' G TZ nj x ■ R\d-u — R\D-v}-> 
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where R spans TZ n> \ and U the set of neighbourhoods of the locus of generic ramification of R. 
This is a fairly coarse topology: for example, the cardinal of the locus of generic ramification 
is not continuous, but only lower semi-continuous in this topology. However, the ramification 
type at O is continuous. On the set X of finite subsets of D*, we put the topology which makes 
the bijection X ~ Un>o((-^*) n ~~ A„)/(5 n a homeomorphism, where A n is the subset of (D*) n 
where at least two components coincide. These topologies do not make the ramification locus 
a continuous function of the ramified covering. Nevertheless, let Ai(1Z n \) denote the space of 
finite Borel measures on 1Z nt \ endowed with the topology of weak convergence. 

Lemma 8.1. The mapping X — > JV[(lZ n \) which sends X to p n \x is continuous. 

Proof — By definition of the topology on X, it suffices to prove that the mapping is continuous 
on (D*) k — Afc for all k > 0. Consider k > 0, X = {x\, . . . ,x k } and a bounded continuous 
function / : TZ n ,\,x — ► ^- Choose e > 0. 

Since lZ n \ x is finite, the continuity of / implies the existence of r > such that the balls 
B(xi,r) are contained in D*, pairwise disjoint for % € {1,. . . , k} and the neighbourhood U = 
B(xi,r) x ... x B(x k ,r) of X in (D*) k - A k satisfies 

VReK nXX ,VR' eO(R,U), \f(R')-f(R)\<e^j . 

Let X' = {x'i, . . . , x' k } be an element of U. Let be a homeomorphism of D such that 4>\d-ij = 
ido-u an d <f>(xi) = x\ for all i G {1, . . . , k}. For each ramified covering tt : R — ► D belonging 
to T^n^Xi the covering &(R) = (4>oir : R — > D) belongs to 1Z n ^ X '- Replacing by its inverse 
in the definition of $ : lZ ni \ )X — > Tl n ,\,x' yields the inverse mapping, hence <E> is a bijection. 
Moreover, the conjugation by (f> determines an isomorphism between Aut(i?) and Aut($(i?)). 
Finally, R and &(R) are isomorphic outside U. Altogether, 

\ P n,xMf) ~ Pn,X,x(f)\ < E #A ut(iZ) lfiHR)) " f(R)l < ^ 

Since k, X, f and e were arbitrary, the result follows. □ 

Let t > be a real number. Let E$ be the distribution of a Poisson point process on D of 
intensity ^ times the Lebesgue measure on D. Under 2$, a random subset of D is contained in 
D* with probability 1 and the average number of points of such a random set is t. Thinking of 
E t as a Borel probability measure on X, we define a measure on lZ n \ by setting 



Pn,A = / Pn,X,X H t (dX). 
JX 



The measure p l nX is finite and satisfies p^ A (l) = e*^) - *. We define a probability measure p} nX 
on TZ ni \ by normalizing p l n x . 

Let us define two functions on 1Z nj \. Firstly, given R £ 7Z nt \, let us define k(R) as the 
number of ramification points of R distinct from O. We have observed that this is a lower 
semi-continuous, hence measurable function of R. Secondly, let x(R) De tb- e Euler characteristic 
of R. 

Lemma 8.2. Let X be a subset of cardinal k of D* . Let R be an element of !Z n ,\,x and 
7 = (<To, • • • , Ofc) a representative of the associated orbit ofV n ,\k- Then 

x(R)=£(a k )-k = £(X)-2d( 7 ). 

In particular, % : 7Z n ^\ — ► Z is upper semi- continuous and measurable. 
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Proof — The first equality follows from the Riemann-Hurwitz formula, the second from the 
definition of the defect of a path. The last assertion follows from the lower semi-continuity of 
k and the fact that £(o~k), which depends only on the ramification type at O, is a continuous 
function of R. □ 



The main result is the following. 

Theorem 8.3. Let N,n > 1 be two integers. Let (Bt)t>o be the Brownian motion on U(N) 
defined in Theorem \3.3[ Let X be a partition of n and a an element of & n which belongs to the 
conjugacy class determined by A. Let t > be a real number. Then 



jnt- 



n 2 t 



We could have avoided the unpleasant exponential factor in the statement of this theorem if 



we had considered the signed measure p t n \ = Jx(~ 1) Pn,A,x ^t(dX) instead of p. 



Proof — Let X be a finite subset of D* . The set TZ nj \,x is m bijection with the set of orbits 
of Hk(o~) under the action of <B n by conjugation. Let R be an element of TZ n ^\ t x and 7 = 
(cr, <7i, . . . , Ofc) a representative of the corresponding orbit. By Lemma [ 



(-l) m PnX x{dR) 



n, A 



Y (_l)M^)^( A )- M (7) 

7en fe ( CT ) 



d>0 



Integrating with respect to 3t(dX), we find 

-l) k (R) N x(R) pl x (dR) 



-ty 



k,d>0 



k\N 2d 



S(a, k, d). 



By Theorem 13.31 the right-hand side of this equality is equal to e * + 2 E 
follows after normalizing p l n \- 



The result 
□ 
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